Gauge theory picture of an ordering transition in a dimer model 
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We study a phase transition in a 3D lattice gauge theory, a "coarse-grained" version of a classical 
dimer model. Duality arguments indicate that the dimer lattice theory should be dual to a XY 
model coupled to a gauge field with geometric frustration. The transition between a Coulomb phase 
with dipolar correlations and a long range ordered columnar phase is understood in terms of a 
Higgs mechanism. Monte Carlo simulations of the dual model indicate a continuous transition with 
exponents close but apparently different from those of the 3d XY model. The continuous nature of 
the transition is confirmed by a flowgram analysis. 
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Realizing that seemingly unrelated phenomena obey 
the same rules is certainly one of the most fascinating 
aspects of physics. In the context of phase transitions 
between different states of matter, this issue of univer- 
sality is often addressed within Landau-Ginzburg- Wilson 
(LGW) theory, where one derives an action in powers of 
the order parameter describing a spontaneous symmetry 
breaking. Recently, new kinds of unconventional phase 
transitions which do not easily fit in the LGW framework 
have been discussed in the domain of strongly correlated 
systems. A first example is the quantum phase tran- 
sition arising in some frustrated antiferromagnets from 
a Neel state with antiferromagnetic order (AF) to a va- 
lence bond solid (VBS) state which breaks lattice symme- 
tries [l[. Another case is the classical interacting dimer 
system on the cubic lattice 0, which displays a contin- 
uous transition from an ordered phase where the dimers 
align in columns to a Coulomb phase with dipolar dimer 
correlations. In all cases, two important issues have to 
been considered: first, one has to make sure that the 
transition is indeed critical and not first-order. Then, 
one must seek for alternative ( "non-LGW" ) descriptions 
of the phase transition. For instance, the possibility of a 
continuous AF-VBS transition has been proposed to be 
understood in terms of spinon deconfinement [l|. How- 
ever, this dcconfincd quantum criticality scenario is con- 
fronted with recent simulations favoring a very weak first- 
order driven process @ . 

In this paper, we focus on the effective description 
of the phase transition in the classical dimer model. 
There, a gauge field arises naturally in the Coulomb 
phase @, |j] . One generally introduces the lattice elec- 
tric field: Ei(r) = (— l) r rij(r) with n being the dimer 
occupation number. To capture the hardcore nature of 
the dimers, a divergence constraint is obeyed: V-E = ±1. 
In the Coulomb phase, the long wavelength properties of 
the system are described by the coarse-grained action: 
S = |r jE 2 d 3 r. Due to the energy form of the micro- 
scopic model, the dimer system is driven at low temper- 
atures to a columnar phase with broken lattice symme- 



tries. Using some duality transformations, we argue in 
this paper that this transition can be understood in terms 
of a Higgs mechanism. At the critical point, the system 
is described by a field theory closely linked to the one ap- 
pearing in studies of deconfined quantum criticality. The 
same scenario has been recently proposed by mapping the 
classical dimer model to a quantum 2d bosonic model 0] . 
We study via Monte Carlo (MC) simulations the transi- 
tion of the effective dual model between Coulomb and 
columnar phases to recover its continuous nature. 

We start by considering a "coarse-grained" model 
where the electric field on the lattice can take all integer 
values, which generalizes the microscopic case. Never- 
theless, we retain the divergence constraint V-E = ±1. 
Our starting point is the following action on the lattice: 
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6 u-j and the angular field 9 is a Lagrange 
multiplier ensuring the divergence constraint. A real- 
valued field is then introduced via the Poisson formula: 
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Performing a duality transformation Q , we introduce the 
dual current qi = Cijk^jUk and the gauge field A defined 
by B = ^CijkVjAk where B. L = E t - Ef. We also de- 
fine the static vector X° by Ef = ^tijk^jXg. The 
circulation of X° is equal to -| modulo 2n on each dual 
plaquette. The dual theory corresponds now to a mag- 
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netic field coupled to the currents q with a static frustra- 
tion field X°. Adding a fugacity term for the currents: 
Sf = jj— J2(li) 2 > one can easu y perform the integration 
over the q variables and find: 

S[Ai, X ] = ^(V^-V,^) 2 

- X m J2^s^7kX-A k -X° k ), (1) 

r* 

with the field x resulting from the conservation of the 
dual currents. The different steps above were already 
discussed for related 3D quantum spin models @, 0| ■ We 
are left with a dual theory of one matter field interacting 
with a non compact gauge field with geometric frustra- 
tion. Note that the non-compactness of the field origi- 
nates from the absence of monomers in the original dimcr 
model. We can now discuss the possible phases encoun- 
tered in this dual model. For small A m , the gauge field is 
essentially free and exhibits dipolar correlations. When 
X, n increases, the matter field condenses and gaps the 
gauge field by the Higgs mechanism. However, this con- 
densation is constrained by the frustration vector X°. 

To find how the matter field condenses, one standard 
possibility [1, 0, 0] is to consider the soft spin version 
of Eq. [1] The kinetic energy turns out to have a two- 
dimensional manifold of minima, giving rise to an effec- 
tive Ginzburg-Landau action function of the two complex 
matter fields 4> — (0i, <p2 ) @: 

S = J d 3 r M(V< - iA t )cf>\ 2 + U(<j>) + y (V x A) 2 ^ , 

(2) 

with C/(0) = r\cf)\ 2 +u|0| 4 +vgls((f)) and where the exact 
form of the 8th order term /§ has been derived in Ref. @ 
for a quantum-mechanical model in 3 + 1 dimensions. 

The criticality of field theories such as Eq. [2] has been 
of great interest in recent years. In (2 + 1) dimensions, 
it is in principle related to the AF-VBS transition at 
the deconfined quantum critical point up to the Is 
term. The connection between the original dimer model 
and the effective action Eq. [2] was recently derived inde- 
pendently ||. Note that the presence of the frustration 
vector X° in Eq. [T] is of crucial importance for the analy- 
sis above. The non-frustrated model, being dual to a XY 
theory, is known to present an inverted XY transition 
which therefore lies in the 3D XY universality class Q. 
Its associated field theory corresponds to a single matter 
field interacting with a gauge field. We therefore expect 
a different behavior with and without frustration. 

Theoretically, the e expansion is of no help to char- 
acterize the transition as it predicts a first order tran- 
sition L 10j for both one component and two component 
matter field model. Several high performance simula- 
tions [ll|, [lH on easy-plane lattice versions of two matter 
fields models have shown evidence for a weak first order 



transition, in contrast with the continuous transition pre- 
dicted by the deconfined criticality scenario. Recent sim- 
ulations on SU(2) versions are controversial, with some 
pointing towards a continuous transition [1 31 ] while so- 
phisticated analyses favor a first-order process Q. 

There are several caveats when considering the two- 
matter fields model. The first problem comes from the 
inherent approximations performed when taking the soft 
spin version of Eq. [TJ and then going back to the lat- 
tice for the numerical simulations. Secondly, simulating 
Eq.[5]is more complicated than Eq.[T]as we have more de- 
grees of freedom with the two-matter fields. Finally and 
most importantly, the connection is lost with the origi- 
nal dimer model, our main motivation. In some cases, 
it might be simpler to go one step back and simulate an 
intermediately-derived field theory such as Eq. [TJ when 
available. This has the advantage of being at the same 
time closer to the microscopic dimer model (and there- 
fore making it possible to use microscopic observables 
such as the columnar order parameter) as well as to a 
gauge-theory description (which is necessary to invoke 
the Higgs mechanism). In this paper, we adopt this strat- 
egy and present results of a MC simulation performed on 
Eq. [1] which shows a transition between a columnar and 
a dipolar phase. The transition is continuous with expo- 
nents found close but possibly different from those of the 
3d XY universality class. We used the Metropolis algo- 
rithm on cubic lattices of size L 3 with periodic boundary 
conditions up to L = 64, taking A m = f3 = 1/T. 

We first take K. = 1 and present thermodynamic re- 
sults. We checked that the probability distribution of 
the action reveals no double peak structure, supporting 
the scenario of a second order transition. The second 
moment of the action M% = ((S — (S)) 2 ) / L 2 displays a 
peak at the transition point (see Fig.QJi). The height of 
the peak grows slowly up to L — 32 and then converges 
within error bars for larger systems. The critical temper- 
ature T^ 2 = 1.05(1) is estimated from the position of the 
maximum. The fact that M% converges at the transition 
indicates a critical exponent a negative, as in the 3d XY 
universality class. We note that M2 appears to converge 
more rapidly with system size than for the non-frustrated 
mod el ( data not shown), similar to the situation found in 
Ref. Il3j. This suggests an exponent a larger (in absolute 
value) than for the 3d XY universality class although it 
is difficult to give a precise value. 

We now consider the low-temperature phase. In this 
lattice gauge theory, a gauge-invariant order parameter 
related to lattice symmetry breaking can be defined. In 
the original microscopic model, dimers order on columns 
as temperature is lowered. Equivalently in the coarse 
grained model, the magnetic field B arranges in staggered 
flux lines in one particular direction. We therefore define 
the local order parameter c;(r) = (— l) r ~ Ti Bi(r) and the 
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FIG. 1: (color online) a) Second moment M2 versus T for 
different system sizes. Inset: Zoom on the peak, b) Columnar 
order parameter C and transverse magnetic field correlator 
G(p min ) as a function of T for different system sizes. 

associated global parameter: 

c = ^lic|| = ^ll£ c ( r )ll- 

r 

At low T, we expect one of the component of the vec- 
tor C to be non-zero, resulting in a finite expecta- 
tion value. The original columnar states of the dimcr 
model are represented by the configurations with: C = 
{±C, 0, 0} , {0, ±C, 0} , {0, 0, ±C}. In Fig. [Tp, C is ob- 
served to vanish at high T, to take a non-zero value as 
T decreases and to finally behave as 1/T at low T as 
can be understood from the equations of motion. To lo- 
cate the critical point, we measure the Binder cumulant 
of the order parameter B = (C 4 )/(C 2 ) 2 , which admits 
a crossing point for different systems sizes (see Fig. [2]). 
This is characteristic of a second order transition and 
leads to an estimate T c co/ = 1.044(5). Assuming the 
standard scaling form B = /(L 1 / iy (T — T c )), the deriva- 
tive dB/dT should scale as L 1 ^ at criticality. We have 
measured this quantity thermodynamically, and display 
its scaling in the left inset of Fig. [2] for the tempera- 
tures T = 1.04 and T = 1.05 around the estimated T c . 
Fits to a power-law form at these two T allow to bound 
0.70 < v < 0.82. The other inset shows the best data 
collapse of the curves according to the scaling form with 
v = 0.73 and T c = 1.044. We find that acceptable data 
collapses can also be obtained for values in the rather 
broad range v = 0.65 — 0.85. 

To characterize the high T phase and detect its 
Coulomb nature, we study the evolution of the trans- 
verse magnetic field correlator at small momentum 

E3: G(p min ) = ^^(PmiJ^iC-Pmin)) with B<(p) 

the Fourier transform of £?j(r) and p min = 2irx/L. 
G(p min ) = 1 if the gauge field presents dipolar corre- 
lations and G(p min ) — ► if the field has short range cor- 
relations. In our case, we expect the field to acquire a 
mass in the low-T phase as the matter field condenses. 
The evolution of G(p min ) with T is shown in Fig. [T] b. 
At high T, G(p min ) — ► 1 and the gauge field is gap- 
less. As T is lowered, G(p min ) decreases and the field 
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FIG. 2: (color online) Binder cumulant of the columnar pa- 
rameter versus T. Left inset: Divergence of the derivative 
dB/dT close to T c versus system size (log-log scale). Right 
inset: Binder cumulant scaling collapse. 
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FIG. 3: (color online) L.G versus T for different system sizes. 
Left inset: Scaling of the derivative LdG/dT versus system 
size L (log-log scale). Right inset: Data collapse of L.G. 

becomes massive due to the Higgs mechanism. Close to 
the critical point, we assume the following scaling ansatz: 
L'w.G(p min ) = /(£ 1/,y (T - T c )), where gauge and scale 
invariance impose t\a = 1 at the critical point [l6j] , The 
crossing of L.G for different system sizes at T c has been 
observed in the non-frustrated model, leading to an esti- 
mate v = 0.67 (l5j |. The frustrated model also presents 
a crossing point (see Fig. [3]) at T C G = 1.035(5). The scal- 
ing versus L of the numerical derivative LdG/dT leads 
to an independent estimate v = 0.73(5) (see left inset 
of Fig. E]), agreeing with the less precise values obtained 
from the Binder cumulant. The value of a obtained with 
hyperscaling a = 2 — vd agrees with the convergence of 
AI-2- The data collapse of L.G obtained with this value 
of v is also of good quality (see right inset). 

So far, if the transition looks well continuous, we can- 
not totally exclude a very weak first-order process, with 
a very large but finite correlation length. In order to 
confirm this, we change the value of the stiffness K, and 
implement the flowgram method 0, [ll[. This method 
relies on the demonstration that the large scale behavior 
for a large value of K, is identical to that at a smaller 
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FIG. 4: (color online) Data collapse of the Flowgram of the 
model. Inset: flows for different values of tC, note how the 
flows look more convergent for smaller stiffness. 



value of JC where the nature of the transition can be eas- 
ily determined. For each JC, we follow Ref. [l3| and define 
the operational critical temperature T C (L) for a size L 
as the temperature where the Binder ratio is equal to 
B c = 1.41 |l7|. We then compute for each value of JC 
and L the value of L.G and recover the flows. If all the 
curves within an interval JC £ [/Ci,/C2] can be collapsed 



into a single master curve by rescaling the system size, 
then it implies that the order of the transition remains 
the same within the interval. A diverging flow for the 
collapse is a clear sign of a first order transition. We 
present here the flowgram obtained for our model with 
JC s [0.36, 2.2] on Fig |4j We have succeeded in perform- 
ing a collapse of the curves by rescaling L — > C{JC)L with 
C(JC) = 2.0//C + 0.4//C 2 . The collapse is clearly converg- 
ing. This suggests that scale invariance is reached and 
goes in favor of a continuous transition in the interval 
considered, confirming the previous analysis at JC = 1. 

To summarize our results, we have shown how, via a 
duality transformation, the interacting dimer model on 
the cubic lattice can be understood in terms of a complex 
matter field coupled to a gauge field with geometric frus- 
tration. As for the microscopic dimer model 0], we found 
a direct continuous transition (by the Higgs mechanism) 
between a dipolar Coulomb phase at low coupling and a 
columnar phase which breaks lattice symmetries at high 
coupling. However, critical exponents are rather differ- 
ent from the ones obtained in the original dimer model 
(where a ~ 0.5 and v ~ 0.5). In particular, the second 
moment of the action clearly shows no diverging peak in 
our case. The anomalous dimension r/A is in agreement 
with the scale dimension of the stiffness in Ref. [H but a 
different value of v is obtained. The exponents that we 
find [y ~ 0.73 and a < 0) seem to bc slightly different 
from those of the 3d XY universality class (and conse- 
quently of the non- frustrated version of Eq. [T]) , although 
we cannot exclude it within our numerical accuracy. In 
both cases, we are left with an interesting open problem 



as there is no field theoretical arguments allowing to say 
that we should end up in the XY universality class. In 
view of our original motivation, an intriguing possibility 
remains that the microscopic dimer model is directly self- 
tuned (for unknown reasons) to a tricritical point. The 
flowgram analysis displays no sign of tricriticality how- 
ever in the range of stiffness considered. Since values of 
JC < 0.36 cannot be reached as the critical temperature 
is too low for the Metropolis algorithm to be efficient, 
one alternative would be to directly perturb the original 
microscopic dimer model to check this scenario. 
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